Derivation of the Matalon-Packter law for Liesegang patterns 
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Theoretical models of the Liesegang phenomena are studied and simple expressions for the spac- 
ing coefficients characterizing the patterns are derived. The emphasis is on displaying the explicite 
dependences on the concentrations of the inner- and the outer-electrolytes. Competing theories (ion- 
product supersaturation, nucleation and droplet growth, induced sol-coagulation) are treated with 
the aim of finding the distinguishing features of the theories. The predictions are compared with 
experiments and the results suggest that the induced sol-coagulation theory is the best candidate 
for describing the experimental observations embodied in the Matalon-Packter law. 



o 

CO 



UGVA DPT 1998/07-1012 



T3 

C 

O 
o 



in 

(N 

i> 
O 

OO 



I 

-a 
c 

o 
o 



X 



I. INTRODUCTION 



Systems that exhibit pattern formation are common in 
nature and patterns often emerge in the wake of a mov- 
ing front [jjj . A well known example that will be studied 
below is the Liesegang phenomenon (||^], the formation 
of precipitation patterns by a moving chemical reaction 
front. In a typical experimental setup, one has a chemi- 
cal reactant (inner electrolyte) dissolved in a gel matrix 
and a second reactant (outer electrolyte) poured onto the 
gel. The concentration of the outer electrolyte is much 
higher than that of the inner one and so it diffuses into 
the gel and reacts with the inner reactant. The reaction 
product then precipitates and, frequently, the observed 
precipitation patterns are a family of bands or rings (de- 
pending on the geometry of the system) clearly separated 
in the direction perpendicular to the diffusion. This pro- 
cess is believed to be responsible for many precipitation 
patterns such as the structure of agate rocks or, to cite 
a somewhat speculative example, the pattern of human 
gallstones. 

Although the Liesegang patterns have been known for 
more than one hundred years ||, there is still disagree- 
ment as to the mechanisms underlying this phenomenon. 
The main problem is the lack of a theory that could 
be compared with experiments quantitatively and, to a 
lesser degree, the lack of experiments that produce re- 
sults amenable to quantitative theoretical analysis. At 
first the task appears to be simple. The Liesegang bands 
are formed at times t n at some positions x n (measured 
from the planes separating the reagents at the initial mo- 
ment), and have a width w n . The quantities t n , x n and 
w n obey the following generic laws: 

i) After a transient time, the position x n of the n-th 



band is related to the time t n of its formation through 
the so-called time law, x n ~ \ft^, first discussed in [0. 
The time law is a rather obvious consequence of the dif- 
fusive motion of a reaction front in a gel and thus it is 
considered to be understood. 

ii) For large n, the ratio of the positions of two consec- 
utive bands approaches a constant value 



l+p 



(1) 



with 0.05 < p < 0.4 typically. This last property is 
known as the Jablczynski law [g| or the spacing law and 
I + p is called the spacing coefficient. Most of the de- 
tailed experimental observations concern this law. It has 
been found that the spacing coefficient is a nonuniversal 
quantity. It depends for example on the concentrations 
ao and bo of the outer and inner electrolytes which can 
be controlled experimentally. Based on several different 
experiments, Matalon and Packter [pP] concluded that: 



P = F(b ) 



G(b ) 
a 



(2) 



where F and G are decreasing functions of their argu- 
ments. For various reagents, they found F(t>o) ~ b 7 
with 0.2 < 7 < 2.7. The function G(bo) is less known but 
it is generally observed that it decreases with increasing 
&o. The observation expressed in equation (^|) is usually 
called the Matalon-Packter law. 

A problem with the Matalon-Packter law is that, ex- 
perimentally, a and b can usually be changed only by 
a factor 5-10 and, furthermore, p itself has a significant 
error bar since it is determined from a finite number of 
x n -s (n ~ 10 — 20). Thus the functional form (||) is far 
from unique, it should be considered as a power law fit 
top 1. 
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iii) Finally, the width law states that the width w n of 
the n-th band is an increasing function of n and typically 
W n ~ x%, a > ||. Since the definition of w n is open 
to debate and the error bars in the measurements of w n 
are rather large, this quantity has largely been ignored 
in the quantitative discussions of Liesegang phenomena. 

Although there are a large number of interesting and 
often puzzling observations about Liesegang phenomena 
in particular systems, the above three points appear to 
be the only ones which describe quantifiable results and 
carry some generality. 

One would expect that so few experimental facts can 
be explained easily. Indeed, there are several competing 
theories [fl0|-[2l| many of which ]To|-[l7|] fare well as far as 
the derivation of the time and spacing laws are concerned 
p2j. All these theories follow how the diffusive reagents 
(ions in the outer and inner electrolytes) A and B turn 
into immobile precipitate D 



A + B 



.c. 



D 



(3) 



by taking into account various scenarios for the inter- 
mediate steps denoted as ...C... . Since the precipitate 
appears through some kind of supersaturation, further 
differences in theories arise from the details of treating 
the nucleation thresholds and the growth kinetics of the 
precipitate. 

The simplest (and first developed) theory is based on 
the concept of supersaturation of ion-product [|10| and has 
been developed by many researchers pl]- p^| , p3| . In this 
theory, there is no intermediate step in-between A, B and 
D. When the local product of the concentrations of the 
reactants, ao, reaches some critical value, q*, nucleation 
of the precipitate D occurs. The nucleated particles grow 
and deplete A and B in their surroundings. As a conse- 
quence, the local level of ab drops and no new nucleation 
takes place. This continues until the reaction zone (where 
ab is maximum) moves far enough that the depletion ef- 
fect of the precipitate becomes weak. Then ab = q* is 
reached again and nucleation can occur. The repetition 
of this process leads to the formation of bands. 

In the next level of complexity, theories contain a single 
intermediate step in ...C... with the mechanism of band 
formation based on the supersaturation of the interme- 
diate compound C [16 17j]. It is assumed that A and B 
react to produce a new species C which also diffuses in 
the gel. The nature of C is not really specified, it may be 
a molecule as well as a colloid particle. When the local 
concentration of C reaches some threshold value, nucle- 
ation occurs and the nucleated particles, D, act as aggre- 
gation seeds. The C particles near to D aggregate to the 
existing droplet (hence become D) provided their local 
concentration is larger than a given aggregation thresh- 
old. These models are characterized by two thresholds, 
one for nucleation and one for droplet growth. The de- 
pletion mechanism is similar to the one described for the 
ion-product theory and it leads to band formation. We 
shall refer to this theory as the theory of nucleation and 
droplet growth. 



A variation of the single intermediate compound the- 
ories is based on the idea of an induced sol- coagulation 
process The compound C is assumed to be the 

sol and this sol coagulates if the following two conditions 
are satisfied. First, the concentration of C exceeds a 
supersaturation threshold c > c* and, second, the local 
concentration of the outer electrolyte is above a threshold 
a > a* that is often referred to as the critical coagulation 
concentration threshold. The band formation is a con- 
sequence of the nucleation and growth of the precipitate 
combined with the motion of the front where a = a*. 

The time- and the spacing laws follow from the above 
theories. Thus, to select the correct theory (provided 
there is a single theory for all Liesegang phenomena) , one 
would have to calculate the functional form of p(ao,6o) 
in order to find agreement or disagreement with the 
Matalon-Packter law. The spacing coefficient, however, 
is obtained from the numerical solution of complicated 
coupled nonlinear differential equations and, as a conse- 
quence, the results are in an implicit form that is not 
particularly useful for deducing p(ao,6o). A notable ex- 
ception is the simplest version of the ion-product the- 
ory which was used by Wagner [filj to derive a result, 

— 2/3 — 

p ~ a , that is at variance with the Matalon-Packter 
law. Wagner's result raises the question whether the 
present theories contain the Matalon-Packter law at all 
and our aim here is to discuss this question and make an 
attempt at answering it. 

In attempting to derive the Matalon-Packter law, our 
basic aim is to keep the calculations and formulas sim- 
ple so that the explicit dependence on ao and bo could 
be displayed. Accordingly, we have to make simplify- 
ing assumptions. The first two of these assumptions are 
about the experimental setup and the majority of the 
experiments satisfy the constraints described there. The 
remaining assumptions are based on observations of the 
time evolution of the concentrations of the reactants and 
reaction products obtained from numerical solution of 
the appropriate equations, from simulations of the pro- 
cess [17j, and also from experience with the analytical 
solution of the problem of chemical reaction zone 
We believe that the important features of the original 
problem are not lost after making the assumptions listed 
below. 

1. The experimental setup can be described by one- 
dimensional reaction-diffusion equations (the gel is 
in a testube whose length is much larger than its 
diameter). 

2. The concentration of the outer electrolyte (ao) is 
assumed to be much larger than that of the inner 
electrolyte (bg) and it is also assumed that a(0, t) is 
kept fixed at x = 0, [a(0,i) — ao]. In experiments 
one typically has 0.005 < bo/ao < 0.1. 

3. Reaction zones and the precipitation bands are 
much narrower than the diffusion length. 
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4. All concentration profiles can be approximated by 
straight lines in the neighborhood of the reaction 
zones and precipitation bands. 

5. The slopes of the straight lines are determined by 
the diffusion lengths and by the asymptotic con- 
centrations of the reacting species. 

6. An existing precipitation band acts as a sink (with 
b(x n ,t > t n ) = or c(x n ,t > t n ) — boundary 
condition) for the inner electrolyte and for the in- 
termediate compound (this is an assumption about 
the reaction and aggregation rates being sufficiently 
large). 

Using the above assumptions, we can derive relatively 
simple expressions for p(a,Q,bo) for all three classes of 
theories (in Sec. 2 for supersturation of ion-product, in 
Sec. 3 for nucleation and droplet growth, and in Sec. 4 for 
induced sol-coagulation theories). The results are com- 
pared both with the corresponding numerical solutions 
of the non-simplified models and with experimental data. 
Sec. 5 contains our final remarks. 

In closing the introduction, we list the symbols fre- 
quently used in the text: 

a(x, t) - concentration of outer electrolyte, A; 
b{x,t) - concentration of inner electrolyte, B; 
a , bo - initial concentrations of A and B; 
k = bo/ao - dimensionless ratio of initial concentra- 
tions; 

c(x, t) - concentration of sol (or intermediate com- 
pound) ; 

D a , D D , D c - diffusion coefficients of A, B, and C; 
Df - effective diffusion coefficient of reaction front 
{A + B^ C); 

q* - ion-product threshold; 

c* - nucleation (coagulation) threshold; 

a* - threshold for induced coagulation; 

x n - position of the n-th precipitation band; 

t n - time of formation of the n-th band; 

1+p = x n+ i/x n - spacing coefficient. 



II. ION-PRODUCT SUPERS ATURATION 

Let us begin with the simplest case of the reaction 
A + B — > D. The ion-product supersaturation theory 
JTo| |l3| is based on the assumption that precipitation of 
D at time t and at a site x occurs if the product of ion 
concentrations a(x,t)b(x,t) exceeds a threshold value 



a(x, t)b(x, t) > q* 



(4) 



^ t = D b^-k6(ab-q*)-Xabd, 
dd 

— = kO{ab - q*) + Xabd, 



(5) 
(6) 
(7) 



where D a and are the diffusion constants of A and B, 
and 9{x) is the step function describing an infinitely sharp 
threshold for precipitation with k being the rate constant. 
The last terms on the right hand sides represent the ag- 
gregation onto the existing precipitate of concentration 
d{x,t). The initial conditions to the above equations 
are given by a(x,0) = a$6(— x), b(x, 0) = bo6(x), and 
d(x,0) = 0. 

These equations can be solved numerically and our an- 
alytical expression found for p(ao, bo) in the above process 
is in accord with the numerical results. Before deriving 
the results, however, let us show how the structure of the 
result may emerge from a simple dimensional analysis. 

If we assume that both the precipitation (fc) and ag- 
gregation (A) rates go to infinity then it follows from di- 
mensional considerations that the dimensionless spacing 
coefficient can be expressed through the available dimen- 
sionless combinations of the parameters and so we have 



p{ b D b q 
a D a a b 



(8) 



Further simplification can be made by assuming that P 
is not singular in its first argument and setting this argu- 
ment to zero since in experiments bo/ao -C 1. Then the 
Matalon-Packter law can emerge only if P is linear in its 
third argument thus giving 



D l Db \ , V I Db \ 1 

D a D a a b 



(9) 



We emphasize that this is not a derivation of the desired 
result. This is, however, what we shall show to be valid 
in the experimentally investigated range (0.05 < p < 0.4) 
of possible p-s. Our calculation is based on a simple an- 
alytical approach combined with numerical evaluation of 
Pi and P 2 . 

Figure [l] gives a characteristic picture showing the con- 
centration profiles as well as the ion-product profile past 
the formation of the precipitation band x n and just be- 
fore the appearance of the (n + l)-th band. It is clear 
from the figure that the linearization of the concentra- 
tion profiles is a good assumption and accordingly, in the 
neighborhood of x n , we can approximate the concentra- 
tions as 



a(x, t) = a>o ^1 
b(x,t) 



b 



v / 2D b (t~t n ) 



{x X n ) 



(10) 
(11) 



In the mean-field theory, the equations describing this 
reaction-diffusion process are given by 



The condition for the n + 1st band to appear is that the 
maximum of the ab product reaches the value q* . This 
leads to the following two equations 
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a(x n+1 ,t n+1 )b(x n+1 ,t n+1 ) = q* 
Q^[a(x,t)b(x,t)] Xn+utn+1 =0 



(12) 
(13) 



The second equation gives a simple relationship between 
t n +i, x n and 



(14) 



\/2D a t n+ i — 2x n +\ — x n 



which allows us to eliminate t n +i from equation (|12j) and 
to obtain an equation for x n +i/x n and x n /x n —±. 



c 0.4 
o 



0.0 



100ab 



100 



150 



200 



250 



FIG. 1. Concentration profiles in the ion-product supersat- 
uration model of Liesegang phenomena. Concentra tions are 
measured in units of ao while the unit of length is */ D a /kao- 
The ionproduct, ab, and its critical value q* are measured in 
units of do and they are magnified by a factor 100 for better 
visibility. 

Assuming now that n is large enough that we can use 
the approximation 



X n +1 
X 1 1 



Xn—X 



= 1+p 



and one obtains the following equation for p: 



Hp) 



? 3 / 2 {i+p) 



(1 + 2p) 2 ^2Tp~ 




(15) 



(16) 



This equation clearly has a p > solution for small val- 
ues of £ = D a q* I (aobo). p increases with £ and 
and goes to infinity for £ — > 1/4, accounting for the fact 
that no bands can be formed if ao&o is too small for ab 
to reach the critical value q*. 

At first sight, the Matalon-Packter law does not fol- 
low from this theory, since p depends on ao and bo only 
through cto&o- Furthermore, another problem is that, for 
small p, we obtain Wagner's result [pT[ , namely the de- 
pendence of p on ao is of the form 



2/3 



(17) 



in contrast to the Matalon-Packter prediction p ~ 1/ao. 
However, one should notice that the range of p where 
p ~ a 2 ^ 3 holds is certainly much smaller than the exper- 
imentally observed values 0.05 < p < 0.4. In this latter 
interval, <&(p) is very well approximated as a straight line 
(seeFig.§) with parameters $(p) = -0.0035+0. 2p. Thus, 
to a good approximation, we can write (16) as a Matalon- 
Packter law with F(b ) = 0.02 and G(V) ~ 1 / b o 



p = 0.02 




(18) 



Pi 



In this approach we obtain P\ = const, and 
^/Db/D a . Numerical solution of the equations (||,^) in- 
dicates that, in reality, the p is given to an excellent ac- 
curacy by 



pw 0.18- 0.052(^) 1/3 

a 



7.5 



aob 
(19) 



For D a — Di,, the order of magnitude for F(bo) ~ 0.12 
agrees with the order of magnitude values of F in experi- 
ments Q and the 1/ao dependence is also the experimen- 
tally observed behavior. However, in experiments, F(bo) 
is not a universal number. It changes not only as Db/D a 
is varied but it also depends on bo [|J. Furthermore, the 
experimentally observed G(bo) is also more varied than 
being just G(bo) ~ l/6o- Thus we must conclude that 
although it is possible to derive the Matalon-Packter law 
from the ion-product supersaturation theory, the result- 
ing functions F(bo) ~ const and G(6 ) ~ l/b are too 
simple to account for experimental observations. 



III. NUCLEATION AND GROWTH 

The theory at the next level of complexity assumes 
that the A + B reaction yields the precipitate D through 
an intermediate compound C. This compound may, in 
principle, be a molecule AB but it may as well be a sol 
particle formed by A and B and, possibly, by some other 
background ions. All this rather complicated situation 
is described by assuming that the C-s can be treated as 
a diffusing species which precipitates if its concentration 
exceeds a threshold c*. The precipitate, D, grows by col- 
lecting the neighboring C-s and the various theories differ 
in the sophistication of the description of this nucleation 
and growth process |l(],[l7]]. In our case the nucleation 
threshold is sharp at c* and precipitate will be assumed 
to be a perfect sink for the C-s. 

In order to understand the formation of precipitate, we 
observe that the reaction zone (where the C-s are pro- 
duced) is narrow |24f| and that c(x, t) reaches its maxi- 
mum there (see FigTpffor characteristic concentration pro- 
files just before the n + 1-st band forms). The n-th band 
acts as a sink and thus the C-s (or about half of them) 
formed at the front, Xf, end up in this band. 
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FIG. 2. Dependence of spacing coefficient 1+p on q*/(ao&o) 
as given by the ion-product supersaturation theory for 
D a = D a = 1. The full line corresponds to the solution of 
equations ([],[]), the dashed line to our approximate theory. 
The dot-dashed line shows the linear dependence for small 
q*/(a b ). 

With a good approximation, we can assume that at 
time i n +i, when c reaches the threshold value c* and the 
(n + l)-th band is about to form, the concentration of 
C-s exhibits a triangular-like shape and it varies linearly 
between c{x n ) — and c(x n +i) = c*. As a consequence, 
we can estimate the current flowing to the n-th band as 



bl = d c 



(20) 



Next we note that the front becomes quasistationary in 
the large-time limit (its velocity goes as i/ ~ l/v^) 
and so we can assume that the above expression for the 
current is valid in the whole time interval [i„,i„ + i] (nu- 
merical solutions of the reaction-diffusion equations sup- 
port this assumption). Then the amount of C which 
disappears into the n-th band during the time interval 
[t n ,t n+ i] can be computed as 



N c pa \j\(t n+1 -t n ) =D c c 



tn-\ 



X n +1 



(21) 



On the other hand, Nc can also be estimated as the 
amount of C-s produced by the front [co(x n +i — x n ), 
where cq is a constant p4j| ] minus the C-s which are in 
the triangle shape, c*{x n +\ — x n ): 

N c » (c - c*) • (x n+1 - x n ) . (22) 

Equating the two estimates of Nc then yields 

D c c* (x n +i ~ x n) 2 



CO 



n+1 



(23) 



We can use now the fact that the band is formed at the 
front whose position, xf, is determined by an effective 
diffusion coefficient Dj j24[] 



Xf(t T: 



n+l 



and find 



l+p/2 



P : 



D f (c ~ c*) 



(24) 



(25) 



The parameters D c and c* in the above expression are 
material parameters. The dependence of p on ao and bo 
enters through the effective diffusion coefficient Dfoi the 
front and through the concentration, cq, of C-s left be- 
hind by the moving front. The remaining task is thus to 
determine Df and Co- We start with the case of D a = Dt 
which can be discussed with relative ease jMj . It has been 
shown for this case that Df/D a is determined from the 
following equation: 



erf 



( D 



\2D a 



1 — K 



(26) 



where erf(x) is the error function |2q| and k = &o/ao- The 
remarkable feature of the solution of this equation is that, 
for the experimentally relevant values 0.005 < n < 0.1, 
the inverse of Df is very well approximated by a linear 
function of k (see Fig.|]) 



Da 

Df 



T]l+T]2-KK 0.158 + 4.03- K 



(27) 



This is an important observation since it follows from 
( p5| ) that p ~ 1/Df for small p and the Matalon-Packter 
law will emerge as a consequence of the above equation. 

The general case D a ^ D is more complicated but the 
results are similar. The effective diffusion coefficient, Df 
is obtained from solution of the following equation pa] 



(28) 



where H{x) = [1 — er f (x)) exp(x 2 ) . The behavior of 
D a /Df as a function of k for D^ = D a = 1 is given on 
Figure ||. We see that D a /Df is linear in k in the experi- 
mentally relevant region and thus equation ( p7| ) emerges 
again but the constants rji and 7y 2 are now dependent on 
D b /D a . 
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FIG. 3. Concentration profiles in the nu- 
cleation-and-growth model of Liesegang phenomena. Units 
are the same as in Fig.l 

We turn now to the calculation of cq. The starting 
point is the observation that the width of the reaction 
zone is much smaller than the diffusion length J24] and 
thus the reaction zone can be approximated as a point 
at Xf(t) = yJIDft where both concentrations a and b 
approach zero. This means that the field b(x, t) satisfies 
the diffusion equation with the reaction term replaced by 
the boundary condition b(x = Xf,t) = and the other 
boundary condition being b(x — ► oo,i) = bo- Once this 
moving boundary problem is solved, the number of C-s, 
Nc, produced up to a given time is obtained by the time 
integral of the current \jb\ — Dbdb/dx evaluated at Xf. 
Finally cq is found by dividing Nc by the advance of the 
front in time t, cq — Nc/xf (here we use the fact that 
the density of C-s produced by the front is constant in 
space |24j). The result of this calculation is given by 



Co 




exp 



El 

2D b 



1-erf 



(29) 



For D a = Db and n -C 1 we can use equation (Pq ) in con- 
junction with the asymptotics of the error function p^ ] 
to obtain cq = b . The physical meaning of this result is 
clear. For n -C 1, we have Db <C Df and thus the front 
moves fast into the region of B-s and, consequently, the 
B-s can be treated as immobile particles yielding cq = bo- 
Corrections to the Co = bo result can also be calculated 
using again the asymptotics of the error function |2(| . To 
first order in Db/Df, we find 



Co 



D f ) °\ D a D f 



(30) 



We can use now the linearity of D a /Df in k — bo/ao to 
obtain the following parametrization for cq: 



Cq = b (<7l + <72K) 



(31) 



where o\ = 1 + r/iDb/ D a and 02 = D a are numbers 

of the order of 1 . Since a± and 01 are of the same order 
of magnitude and since k <C 1, in the following we shall 
neglect the k term in (pHl) 




FIG. 4. Inverse of the effective diffusion constant of the 
front Df as a function of k = bo/ao, for D a = Db = 1. 

It is important to note here that a similar omission of 
the bo/a-o term would not be justified in the linearized 
form of D a /D f @. There we have a constant term, 
■qi 0.2 that is an order of magnitude smaller than co- 
efficient, 772 ~ 4, in front of bo/ao- Thus, for the relevant 
values of bo/ao, the two terms contribute equally. 

The expressions ( pti| , pT| , p2[ ) remain good approxima- 
tions up to the point where Db/Df « 1. This happens, 
however, only at rather large values of D b /D a for n <C 1. 
Since the diffusion coefficients of usual electrolytes in 
aqueous solutions are usually within a factor 2 of each 
other |27j and their ratios rarely exceed 5 |^8| , we believe 
that the approximation ( j3ll) can be used for realistic ex- 
perimental setups. 



o Simulations of the R-D equations 
— Theory with c =c (b ,K) 
-— Theory with c o =c o (b ,K=0) 
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FIG. 5. Spacing coefficient for the nucleation-and-growth 
theory. The k = bo/ao dependence of p is shown for the 
following choice of parameters: Db/D a = 1, D c /D a = 0.2, 
c* /bo — 0.633 and with all the reaction rates taken to be 
large. The dashed line is the "linear" Matalon-Packter law 
( ^3[ ) while the solid line is the "non-linear" version of Jjil ) 
where the K-dependence of Co (0) is kept. 



c « ax b 



(32) 



Having determined Df and cq, we can now return 



G 



to equation fl25| ) for p and find again a version of the 
Matalon-Packter law: 



D C C*T)\ 



+ 



D c c*r] 2 b 



D a (o"i b - c* ) D a (ax b - c* )o 



G(bp) 
a 

(33) 



In order to gain confidence in the above result, we have 
calculated p for k in the experimental range 0.005 < k < 
0.1 by solving using the appropriate reaction-diffusion 
equations [equations (§|||@) must be modified and sup- 
plemented by another equation in order to take into ac- 
count the nucleation and growth processes described at 
the start of the section (for details see Jl7j|)]. The results 
are displayed on Figure || where the Matalon-Packter law 
derived above (solid line on the figure) is shown to per- 
form very well considering the simplicity of the deriva- 
tion. It should be noted that nonlinear dependence on k 
sets in for k > 0.06 and that taking into account the k- 
dependence of cq (dashed line) modifies the straight line 
with the right curvature. 

One can observe from equation ( |33] ) that the 
nucleation-and-growth theory gives more complicated bo 
dependences for F(bo) and G(&o) than the ion- production 
theory. Power law forms for F and G are found in the 
limit of c* <C bo (which might be the experimental limit) 
where one obtains F(bo) ~ l/£>o an d G(bo) ~ constant. 
It is remarkable that these are similar to the ion-product 
results but with the roles of F and G interchanged. Com- 
paring the two theories, we see that the nucleation-and- 
growth theory performs better in the sense that it can 
produce power law behavior for F(bo) as observed in 
some experiments and, at the same time, has a G(bo) 
which depends on 6 more weakly but it is a decreasing 
function of bo in agreement with experiments. 

In conclusion, the nucleation-and-growth theory pro- 
vides us with a Matalon-Packter law that is closer to 
experiments and, perhaps, describes some of them. Ne- 
vertheless, the functions F and G appear to be too rigid 
to accomodate all experimental findings. It should also 
be noted that this theory has a prediction G(bo) / F(bo) ~ 
bo that is experimentally easily distinguishable from the 
prediction of the ion-product theory G(bo)/F(bo) ~ 

l/bo. 



IV. INDUCED SOL-COAGULATION 



The induced sol- coagulation theory |14| , |15| | is a general- 
ization of the nucleation-and-growth theory The sol, C, 
is produced at the reaction front in the A + B — > C reac- 
tion and it flocculates if the following two conditions are 
satisfied. First, c must exceed a supersaturation thresh- 
old c* and, second, the concentration of the outer elec- 
trolyte, a, must be above the critical coagulation con- 
centration threshold, a*. The second condition arises in 
systems where the A ions screen the repulsive electro- 
static interaction among the sol particles. 



Fig.|6| shows characteristic concentration profiles in this 
process just as the n + 1-st band is appearing. The 
first remarkable feature of this picture is that the re- 
action front is way ahead of the precipitation zone (this 
does seem to happen in some experiments Jl9|| ). This 
is understandable, since a(xf,t) — > at the front while 
at the place where precipitation occurs one must have 
a(x n+ i, t n +i) > a* . The second important feature is that 
the reaction front leaves C-s behind at a fixed concen- 
tration Co, as already discussed in the previous section. 
In order that precipitation could occur, the parameters 
must be assumed such that Co > c*. Further assuming 
that the sol does not diffuse too fast (D c -C D a ), one of 
the condition for precipitation (c > c* ) is always satisfied 
behind the front and far away from the last existing band 
(x„). 

Consequently, the position, x n +i, where the next band 
appears will be determined by the arrival of the con- 
centration 'front' a = a* to a position where c = c*. 
Since the a and c profiles near x n can be written as 
a = ao(l — x/xf) and c = co(x — x n )/-^2D c (t — t n ), 
the above conditions yield 



Cp{x n+ i - X n ) 

y/2D c (t n +i — t n ) 

Xf(t n+ l) 



a 1 



(34) 
(35) 



where Xf (t) = ■J2Dji is the position of the reaction front 
at time t. 

In order to calculate p, we note now that equation ( |35| ) 
yields 



Xf(tn-\ 



«() 



Xf(t n ) 



(36) 



The above equation can then be used in conjunction with 
Xf(t n +i) = y/2Dftn to obtain p from equation ( p4] ) 



l+p/2 



P 



2D c c* 



D f 4(l-a*/a y 



(37) 



Since the values of Df and cq come from reaction A+B — > 
C which is decoupled from the the next stages of nu- 
cleation and growth processes, the calculations and ap- 
proximations of these quantities can be taken over from 
Section III. Moreover, we shall assume that a* is much 
smaller than ao so that the right-hand side of eq. (^) can 
be expanded in a* /ao- There is no experimental informa- 
tion on a* and, although a* /ao -C 1 seems to be a natu- 
ral assumption, we must be cautious with the a* /ao — ► 
limit. The problem is that the picture we are working 
with for the induced sol-coagulation theory (Figure |^) is 
not valid for arbitrary small a* . Indeed, the fact that the 
position of the A + B — > C reaction front, Xf is far ahead 
of the point of formation of the n + 1-st band means 
that diffusion length of the C particles is smaller than 
Xf(t n+ i) - x n+ i 
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y/2D c {x n+1 - X n ) < Xf(t 

n+1 ) 



(38) 




300 



FIG. 6. Concentration profiles and thresholds in the in- 
duced sol-coagulation theory. Units are the same as in Fig.l. 



Using equations Q3gj37|) and neglecting higher order 
terms in p and a* /an, the above equation yields 



D c c* ap 1 
D f c Q a* 2 



(39) 



The meaning of this result is the following. The place 
where the new band nucleates moves to the A + B — > C 
reaction zone and thus the inequality ( |38| ) gets violated 
upon increasing D c or c* and decreasing a*. It is clear 
thus that, at fixed values of the other parameters, a* /a® 
cannot be taken to zero. 



that the result for the nucleation-and-growth theory ( J35 
cannot be obtained from (|37|). The matter of finding an 
expression for p valid in the a* / 'an — > limit and reducing 
to ( ^3|) appears to be a rather difficult task since various 
lengthscales become comparable approaching this limit. 
We have not been able to obtain such an interpolating 
formula. 

Returning to (|39|), we see that the material parame- 
ters can be such (e.g. if c*/co <C 1) that a*/ao can be 
small without the inequality ([59]) being violated. For 
such range of parameters we can expand (p7|) in a*/ao 
and obtain again a version of the Matalon-Packter law: 



2D c c* 2 r h , 2C c c* 2 (2 m a* + V2 b Q ) 



D a a\ b 2 D a a\ 



hi 



'5 a o 



(40) 



The spacing coefficient ( f!o| ) can again be compared 
with that obtained from the appropriate reaction dif- 
fusion equations (the equations for the nucleation-and- 
growth theory must be augmented with the condition 
that nucleation can occur only if a > a*). An exhaus- 
tive numerical study is practically impossible due to the 
number of parameters in the problem (k, a* /an, c*/ao, 
Db/D a , D c /D a ), plus the rates of reaction, nucleation, 
and aggregation). One can easily miss regimes of non- 
trivial behavior in this high dimensional parameter space 
and we claim with this numerical study only that the 
Matalon-Packter law as given by equation ( [io| ) is indeed 
observed for a reasonable range of parameters (Figures 

As one can see from (^OJ), the induced sol-coagulation 
theory provides us with 



Theory with c =c (b ,K=0) 
Theory with c =c (b ,K) 
Simulations of the R-D equations 



1.35 



+ 1.25 



I 



I 



I 



1.15 - I 



1.05 



0.00 



0.04 



0.08 



FIG. 7. Spacing coefficient for the induced soi-coagulation 
theory. The dependence of p on k, — bo /an is displayed for the 
following choice of parameters: Db/D a = 1, D c /D a = 0.3, 
c* /bo — 0.86, a* /bo = 7.8, and with all the reaction rates 
taken to be large. The dashed line is the "linear" Mat- 
alon-Packter law (^) while the solid line is the "non-linear" 
version of ( fib] ) where the K-dependence of Co (^l|) is kept when 
going from eq.(l37|) to eq.(EoT). 



The lack of the limit a* /an — > has the consequence 



F(b Q )~l/b 2 Q , G(& )~(a + A))A>o 



(41) 



The power law form F{po) 



is close to what 



has been observed in some experiments ]6j and the fact 
that G(6o) is a decreasing function of bo is also in ac- 
cord with the observations. It should also be clear that 
a more precise theory of the induced sol-coagulation pro- 
cess should reproduce, in the a*/ao — > limit, the re- 
sult F ~ b^ 1 obtained in the nucleation-and-growth the- 
ory. Thus, assuming that the crossover between the b^ 2 
and the b^ 1 behaviours is smooth, one should be able 
to find regimes where F(bo) ~ 6q with 1 < 7 < 2. Since 
this covers a large portion of the experimentally observed 
range 0.2 < 7 < 2.7, we conclude that the induced sol- 
coagulation theory provides the best description of the 
pattern formation in Liesegang phenomena. 



V. FINAL REMARKS 

Our main results are summarized in equations (|lS|), 
(|33|), and ( f40| ) giving the Matalon-Packter law for the 
three main theories, respectively. Comparing these for- 
mulas, we arrived at our main conclusion, namely that 
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the induced sol-coagulation theory is the best in describ- 
ing the experimental observations on the spacing coeffi- 
cient of Liesegang patterns. 

It is no doubt that the arguments used in the deriva- 
tion of p can be refined and made more precise. The main 
aim of our work, however, was the demonstration that the 
Matalon-Packter law can be understood in terms of sim- 
ple pictures for which it is possible to develop analytical 
arguments. We hope that one can build on these results 
and achieve a better understanding of Liesegang phe- 
nomena. In particular, it would be important to find a 
description of the induced sol-coagulation process which 
contained the crossover to results for the nucleation-and- 
growth model. Furthermore, we have just considered the 
simplest cases associated with the reaction A + B — ► C . 
In experiments, equally important are the A<i + B — > C 
and A + B2 — > C cases and, obviously, one should also 
explore the possibilites associated with these and, pos- 
sibly, more complicated reaction schemes. Finally, an 
interesting and important test of the conclusion that the 
induced sol-coagulation theory is preferable in the studies 
of Liesegang phenomena would be the application of this 
theory for the quantitative description of revert patterns 
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